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Challenges for MAV 


• Low Re (100 to 10 5 ) 
aerodynamics. 

- Three-dimensional 
geometry 

- Unsteady flow 

- Laminar separations 

- Instabilities 


Re ~ 100 

AeroVironment 
Black Widow MAV 


Re ~ 100,000 



Re ~ 10,000 
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• Closed-loop control 

- Nonlinearity 

- High-dimensional 
systems. 


- Feedback to achieve 
flow states not available 
in open loop 
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Opportunities 



Birds, insects, and bats offer enticing performance 
benchmarks: 

- hover, low cruising speed 

- agility 

- propulsive efficiency 

- robustness, gust resistance 

- sensor-based control 
High-lift: stable leading edge vortex 
Maneuverability: control/synchronization of vortex 
formation/shedding 



Can closed-loop control be used to 
mimic some of these behaviors on a 
(more) conventional airfoil? 





00 
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MURI 


Simulation and 
Reduced-order 
modeling 


Unsteady wind tunnel 
experiments at 
moderate Re 



Oil tunnel experiments 
at low Re 




Control law development: 

> Stabilize LEV (enhance lift) 

> Suppress vortex shedding 

> Synchronize vortex 
shedding 

> Lift cancellation in gusting 
flow 

> Maneuver (perched 
landing) 
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LEV control, semi-circular wing, Re=68000 



Pulsed-blowing actuation concentrates 
vorticity at leading edge. 

F + = fc/U = 1.1 
C M = .0074 
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LEV Control - Rectangular AR=2 wing, Re=300 
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Numerical Simulation (DNS) 


• Re 0(1 0 2 ) to 0(1 0 3 ) 

• Data for modeling efforts & test 
control in simulation 

• Validated Immersed boundary 
projection method 

• Fast algorithm (FFT+multi- 
domain far-field BC) 

• Wrappers for: 

• Linearized/adjoint 

• Linear stability 

• Continuation (stable/unstable 
steady states) 

• Snapshots for 
POD/BPOD/Galerkin 
Projection 
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Model problem: 2D flat plate at angle of attack 


• Economical test case 
displaying vortex 
shedding behavior 

• Close to bifurcation point 
(with Re or AOA) 

• Use actuation and 
feedback to: 

- Stabilized otherwise 
unstable equilibria 

- Phase-lock lift 
fluctuations to actuation 

- Actuation modeled as 
localized body force 
near leading or trailing 
edge 



Leading edge actuator 



Trailing edge actuator 


* 8 


00 
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The flow becomes phase-locked to the actuation upto 15deg, but at higher alpha, it displays a more 
complicated limit cycle, at which subharmonic to the forcing frequency is excited. 

PINK stars represent the period-averaged lift over each actuation period. 

And they have the same value if the flow is phase-locked, as in the case of lower alpha. 


Lift history over time on the RHS show the example of these subharmonic behaviors at alpha=50 

Each subharmonic limit cycle is consisted of several actuation periods associated with a particular 
phase shift between the forcing signal and the lift, and with different period-averaged values, as 
shown as variation of PINK stars on the LHS for high alpha. 
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Phase lock loop 

• Can we obtain phase lock of ‘good’ limit cycles without 
searching for the right forcing frequency? 

• Can we achieve lock on starting from arbitrary 1C (or in 
presence of noise)? 

• Adjust (slightly) forcing frequency and phase to attempt 
to match specified phase shift with output lift signal 


Uj et (t) = a 0 + K p [ai(t) cos(cjit + <ft) + bi(t) sin^j t + fa)] 



Feedback Signal 



SO 

o 
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Result 



o r ■ : 1 « -L- 

10 20 30 40 50 

AOA 


Base 

Open-loop, co f = 1 

Feedback 


• Achieved high-lift unsteady flow states that 
open-loop control can neither achieve nor 
sustain 


* 


In order to phase lock the flow at the desired shedding cycle, particularly at Phi, best, 
We designed a feedback compensator. 


(Even though the open-loop forcing at Wf below Wn can lead to phase-locked limit cycles with a high 
average lift,) 


This feedback controller resulted in the phase-locked limit cycles that the open-loop control could not 
achieve for alpha=30 and 40 


Particularly for alpha=40, the feedback was able to stabilize the limit cycle that was not stable with 
any of the open-loop periodic forcing. 

This results in stable phase-locked limit cycles for a larger range of forcing frequencies than the open- 
loop control. 

Also, it was shown that the feedback achieved the high-lift unsteady flow states that open-loop control 
could not sustain even after the states have been achieved for a long period of time. 
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Optimized waveform 

• What is the ideal control (<|> = U jet ) to maximize lift? 



Determines “direction” to 
change control to 
decrease J 


Determines magnitude 
of change 


j = i 


1 , 


to 


f 2 (A(t); x; t) dxdt + 

body 

Maximize lift 


A 2 (t)dxdt 


t 0 body 

Penalty to boun4 i 
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Receding horizon approach 




• Problems get exponentially 
harder to optimize as 
prediction horizon is 
increased. 

• Controls near the end of each 
optimization horizon 
discarded. 

• Chess: re-evaluation of the 
game plan after each move 
played. 


Specific example 
achieving (about) the 
same period-averaged lift 
as sinusoidal, but with 
40% of the energy 
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Optimal control is still essentially open loop 

• Apply Optimized Waveform in Open loop to Baseline 



Time 




Optimization provides a periodic control waveform that achieves period limit cycle after 4~5 transient 
periods. 

However, it’d be practically implementable only if we can reproduce the high lift limit cycles without 
the transient signal started from any phase of the baseline cycle. 


While it is straight forward to extract a single period of the optimal waveform 


If optimized waveform is applied in open loop to the baseline flow, the flow fails to lock to the forcing 
and the performance can be significantly degraded as shown in the figure. 


This calls for designing a feedback algorithm 

To achieve phase lock between the lift and the optimal waveform. 


14 
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We feedback lift again as an attempt to march along optimized Ujet accordingly. 

The goal is to shift and deform optimal Ujet with consistent phase difference between each 
harmonics. 

However, this is a complex task since optimal Ujet consists of more harmonics than a pure sinusoidal. 
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Recall from the previous implementation of feedback with sinusoidal waveform; 

feedback controller demodulated the lift signal, applied a low-pass filter, shifted the phase by a 
specific amount, and remodulated the signal in order to produce a sinusoidal output locked by a 
specified lag to the lift signal. 


This can be used as a narrowband filter if no phase shift is added. 

The output signal, y(t) is now a sinusoidal which retains the dominant frequency of CL , coo and filters 
out higher harmonics. 

Since we’re interested in marching optimal Ujet at the dominant frequency of CL , it is easier to track 
the frequency of y(t) instead of CL . 
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we use Extended Kalman Filter (EKF) to dynamically estimate frequency Wo, and the phase theta(t). 
Based on the phase estimate theta(t), we can then march along optimal Ujet as shown. 

Where theta, desired is an additional (specific) phase shift relative to the lift signal, (theta, desired=0, in 
optimization case) 
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Also when computing y(t), note that if demodulation frequency, coi is not equal to the actual frequency 
of the lift signal, coo , then y(t) will not be in phase with CL (t), 

But they will be in phase only when ooi = coo . 

Thus, it is necessary to update coi periodically with “ coo , the frequency estimated by an EKF. 

Now, y(t) will be in phase with CL (t) with the same dominant frequency. 

Based on the estimated phase ~ 0, we can then march along the optimized waveform as, 

To summarize : 

1 . Narrowband filter is used on the lift cycle to obtain a sinusoidal signal. 

2. Filtered lift signal is used as input to frequency tracking Extended Kalman Filter (EKF) to estimate 
phase, ~ 0(t) of the lift signal. 

3. EKF frequency estimate is used to tune the filter to avoid introducing phase lag. 

4. Finally, phase estimate * 0(t) from EKF is used to march along cpoptimal . 
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Feedback Control w / Optimal U jet 



• Feedback achieves the same limit cycle as the optimization 
started from any phase of the baseline shedding cycle 


• Feedback allows us to phase lock an arbitrary waveform 
at any desired phase shift. 


• Can study which features are critical to high lift. 

^Tit- 


Now we can start from any phase of the baseline flow, 

and the feedback (in red) phase-locks the flow at the same phase shift, and achieves the same high- 
lift limit cycle as the optimization. 


The feedback controller now allows us to phase-lock an essentially arbitrary waveform at any desired 
phase shift, 

Thus we can utilize this fact to investigate which features of the optimized waveform are critical to 
high lift. 



NASA/CP — 20 10-216112 200 


Stabilize unstable equlibrium (optimal control)? 



^ = f [ (f y (<j>(t),x,t)) 2 dxdt + C w [ [ (4>(t)) 2 dxdt 

J to J Q “to J 


Minimize lift 


□ * 
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Stabilize Steady State at a = 15deg 


LEact. MinRMS, aoa=1 5deg, [Cw=.01 , Fmax = 0.0, Ci = 0.0] 
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Model-based feedback control using low-order models 

S. Ahuja and C. Rowley, Princeton University 


10 mode BPOD model 
Full state feedback (also 
possible with 2 
observers) 

■ LQR to determine K 

■ Large domain of 
attraction even in the 
full NL system 
Controller suppresses 
the vortex shedding 
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Summary 


• Open loop LEV control shows lift enhancement and flow 
structure similar to dynamic stall vortex 

- Extra lift used with dynamic physics-based models to cancel lift 
fluctuations in gusting flow (D. Williams) 

• 2D vortex shedding 

- Post-stall open-loop forcing gives complex, subharmonic resonance 

- Individual periods of forcing appear to be unstable periodic orbits 

- Phase Lock Loop able to phase lock lift fluctuations to forcing, 
stabilize periodic orbits (e.g. max lift) 

- Optimization of waveforms using gradient-based optimal control, 
implemented in closed-loop with PLL 

- Supression of vortex shedding with BPOD/linear 
instability/LQR/observer models 

• Extension of the BPOD approach to unstable period orbits 



23 
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Model-based feedback control using low-order 
models 

• Computationally tractable implementation of balanced 
truncation 

• Balanced POD based on snapshots of impulse (forward) 
response, adjoint response 

• Bi-orthogonal set of forward/adjoint modes 

• Galerkin Projection, retain small number of modes 

• A priori error bounds for the ROM 

• Extension to unstable equilibria 

- BPOD only on stable manifold 

- Unstable manifold (2 modes) treated exactly 

- Similar extension to periodic orbits 


